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Abstract 

We present a theoretical discovery of an unconventional mechanism of inverse Faraday effect 
(IFE) which acts selectively on topological magnetic structures. The effect, topological inverse 
Faraday effect (TIFF), is induced by spin Berry's phase of the magnetic structure when a circularly 
polarized light is applied. Thus a spin-orbit interaction is not necessary unlike in the conventional 
IFE. We demonstrate by numerical simulation that TIFF realizes ultrafast switching of a magnetic 
vortex within a switching time of 150 ps without magnetic field. 



1 



Ultrafast magnetization switching without magnetic field is required from the view point 
of spintronics apphcation for achieving high speed read- write processes in non- volatile memo- 
ries. Current-induced magnetization dynamics is one of the possible phenomena for realizing 
fast processes. By applying a charge current to the magnetic domain wall, the domain wall 
moves with the velocity of ~ 100 m/s-, which corresponds to a switching time of 1 ns when 
the device size is 100 nm. A magnetic vortex core is another candidate of a memory device 
where the polarization direction of the core is used as information. If one uses a magnetic 
field, switching of a vortex core can be carried out within 140 ps -, but use of magnetic 
field has a disadvantage in realizing high density memories. A switching by field-driven spin 
wave excitation was recently carried out and switching time of 200 ps was obtained^. In a 
switching of a vortex core carried out by applying an electric current pulse, the switching 
time was ~ 20 ns There are, however, two difficulties to apply the current-induced vortex 
core switching to memory devices. First, a high current density (~10^^ A/m^) is needed. 
Heating effect due to high current density causes an unstable operation. Second, a slow 
processional motion of the core is needed before the core switching occurs. This transient 
motion makes it difficult to predict when the core switching occurs. It is therefore desired 
to find other mechanisms for realizing fast vortex core switching. 

Recent studies using femtosecond laser pulses revealed fast magnetic switchings^"-. A 
subpicosecond demagnetization of Ni was induced by a heating effect of a femtosecond laser 
pulse^. In Ref. 6, a magnetic field was applied besides a laser pulse to realize magnetiza- 
tion reversal in TbFeCo. Recently, thermal reversal of a ferrimagnet GdFeCo by applying 
a laser pulse was carried out, and reversal mechanism was explored^. Use of circularly 
polarized light^^i^ has a great advantage in realizing magnetization switching in terms of 
all-optical methods without magnetic field. Circularly polarized light induces a magnetiza- 
tion switching by inverse Faraday effect (IFE)^'^^. An effective magnetic field generated by 
the polarized light via conventional IFE turned out to reach 20 T for a strong laser pulse, 
and such effective field makes possible a switching in a time scale of 60 ps -^2.. The con- 
ventional IFE emerges in the presence of uniform magnetization and spin-orbit interaction. 
To decrease the switching time, therefore, one needs to use heavy elements to enhance the 
spin-orbit interaction. 

Recent studies on the anomalous Hall effect have pointed out the similarity of spin-orbit 
interaction and topological magnetic structuresi^*^. In fact, topological structures carrying 



2 



spin Berry's phase couple the coordinate space with spin space just hke spin-orbit interaction 
does. One can thus expect that the role of the spin-orbit interaction in IFE can be replaced 
by spin Berry's phase of topological magnetic structures. A great advantage of using the 
topological magnetic structure is that the effective flux of the Berry's phase can be as large 
as one flux quanta = 4 x 10~^^ Wb) per unit area of lattice. For a vortex with core radius 
of 10 nm, the effective field of spin Berry's phase thus reaches 6 T. 

In this letter, we show that there is indeed an unconventional IFE acting selectively on 
topological magnetic structures with spin Berry's phase or chirality, such as magnetic vortices 
and skyrmionsi^. By analytical calculation based on the nonequilibrium Green's functions, 
we derive a general formula of the effective magnetic field induced by the incident polarized 
light and the Berry's phase of magnetization structure. We call the present Berry's phase 
induced IFE as topological IFE (TIFE). The effect is useful to control individual vortex 
or skyrmion without touching the background, and has advantage for designing devices 
based on topological magnetic structures. Furthermore, we demonstrate that the ultrafast 
magnetic vortex core switching in 150 ps can be achieved in magnetic disk by using a 
numerical simulation. 

The problem we are going to solve is to construct a vector representing the effective 
magnetic field induced by TIFE, i/xiFE, from three vectors, L proportional to the helicity 
of incident light, direction of local spin n and Berry's phase field, From a symmetry 
argument, we have three possibilities; i/xiFE oc n ■ X), X ($ ■ n), and $ (L ■ n). Which 
coupling emerges in reality is answered by a microscopic calculation. We will show that 
-f^TiFE oc n ($ ■ 2y) is the correct answer. The effective field indicates that the energy density 
of the system when a circularly polarized light is injected is -Etife = — fi'/^B^'-ffTiFE ■ oc 
($ ■ L), where g is the Lande factor, /^b is the Bohr magneton, and S is the magnitude of 
local spin. The angular momentum of light is therefore coupled to the effective magnetic 
flux of the spin Berry's phase, in the same manner as the orbital (£) or spin {S) angular 
momentum couples to a magnetic field {H) as I ■ H and S ■ H. 

Let us derive an analytical expression for TIFE by a microscopic calculation. We consider 
conduction electrons interacting with a local spin structure by a s-d type interaction. The 
electric field of the incident light is included to the second order to describe the circular 
polarization. We consider a light in the THz regime. This simplifies theoretical study greatly, 
since then the electron excitation is limited only very close to the Fermi level. Nevertheless, 
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the results obtained in the present paper are expected to be quahtatively extended to the 
case of the visible light, as suggested by qualitatively the same results for visible light and 
THz light in the case of conventional IFEj^^. 

The spin structure is represented by static classical vectors, S{r) = Sn{r), where n 
is the unit vector representing the direction. The s-d interaction between the conduction 
electron and localized spins reads 

-Hsd = -^j d'^rn ■ (cVc) , (1) 

where c = {c^, ci) and are the annihilation and creation operators of the conduction 
electrons, respectively (indices ^ and 1- represent spin). We consider a strong coupling case, 
where the s-d exchange splitting. A, is large. The conduction electron spin thus follows 
almost perfectly the localized spin structure, namely it is in the adiabatic regime^^. The 
electric field of applied light is represented by -E = Re ^Se^^^''^~^^^ , where £ is complex 
amplitude, and Q and Q are wave vector and frequency, respectively. 

We use a spin gauge transformation to extract low energy excitations.— Defining an elec- 
tron operator a in the gauge transformed frame as a(r) = U{r)c{r), where ?7 is a 2 x 2 
unitary matrix, we diagonalize the s-d interaction as Hsd = J d'^ra^aza. This diago- 
nalization is accomplished by choosing U = m ■ cr with m = (sin | cos 0, sin | sin 0, cos |) , 
where 9 and <p are the polar coordinates of n —. By the above local gauge transformation, 
an SU(2) gauge field, Ag • = — |Tr [a°'WdiU~\ (a = x,y,z), arises which couples to the spin 
current, defined as jg ^ = ■^[{dfj_a^)a°'a — a^a'^d^a — iA'^^a^a]. The total Lagrangian for the 
conduction electron thus reads 

d 



C = I d'ra^ { ih— + — + ep - Aa^ 
ot 2m 
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where ji = -^[o-^dia — {dia))a — ij^AiO^a + 2iy4g j(aV"a)] is a charge current, which couples 
to a vector potential of light. A, defined as -E = —dtA. 

The TIFE is studied by calculating the local spin density induced by the incident light, 
s{r) = {c^crc). In terms of the lesser Green's function^^ for the gauge-transformed opera- 
tor, G'^{r,t,r',t') = j^{a){r\t')a{r,t)), z-component {i = x,y,z) of the spin density reads 
Sj(r) = —ih'^j7lijTT[a^G'^{r,t,r,t)], where TZij = 2mimj — Sij represents components of 
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FIG. 1: (color online) (a) Feynman diagrams representing the effective magnetic field generated 
by an topological inverse Faraday effect (TIFF). Solid line represents electron's Green's function, 
and wavy line represents the electric field {£) of the incident circularly polarized light. Contribu- 
tions (al) and (a2) are linear and second order in spin gauge field (As, denoted by dotted line), 
respectively, (b) Schematic depiction of the TIFF induced by a circularly polarization, £ x f *, 
and spin chirality, The induced magnetic field, i/xiFEi is parallel or anti-parallel to the local 
spin direction, n. (c) Profile of effective magnetic field, i^TiFE (Eq. dSl)), induced for a vortex as 
a function of distance from the core center, r. The light intensity is L = 0.3 x 10~^® Tm^. 

a 3 X 3 rotation matrix. The effective magnetic field induced by TIFE thus 



The lesser Green's function is evaluated by a standard perturbative method^. The 
lowest-order contributions to TIFE are diagrammatically shown in Fig. 1 (a). The vector 
potential of the incident light, A, and the spin gauge field. As, are included to the second 
order. The contributions are calculated expanding Green's functions with respect to the 
frequency of THz light, fi, and to the wave vector of the spin gauge field, q, and we keep 
the lowest contributions. We note here that the contributions linear in As (the first two 
diagrams of Fig. 1 (a)) and the one second order in As (the last four diagrams of Fig. 1 
(a)) result in the contribution of the same form proportional to Cj^m^^^sm owing to the 
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identity satisfied by the spin gauge field^^, deA^^^ — dmA^i — ^iS-yS^I ^A^^^ = 0. After a 
straiglit forward calculation, we obtain the lowest contribution to the effective field as 

i^TIFE, = 5^«^,(^^Xr)., (4) 
3 

where = iaY^em'^iz^jemdeAl^,,, and a = (i^t + (^t + uf (^F^rf - eptr-^^)- 

Here ^'f-(^) is the density of states, t^(^\^) is the electron's lifetime, and eFt(i) is the Fermi 
energy, which depend on the spins (up or down), respectively. Coefficient a^j is simplified 
using spin chirality field defined as $ j = 'Yltm ^iemdeA^^ and by noting TZiz = Ui asr^ 

aij = iaui^j. (5) 

In terms of vector n, scalar chirality of a spin texture is $j = ^^-^ €ijkn ■ {djU x dkn), where 
the direction of $ is perpendicular to the plane the spins lie. Equation (jS]) indicates that an 
effective field arises when the magnetic structure has a finite spin chirality ($) like vortices 
and skyrimons. The effective field is also written as 

J/tife = {L-^) n, (6) 

where 

L = m{£x £*) , (7) 

is proportional to the light intensity and helicity. As seen in Eq. ([6]), the direction of the 
field is parallel or anti-parrallel to the local spin direction, n (see Fig. 1(b)). The magnitude 
of the effective field is determined by the angle between the helicity of the light, i{£ x £*), 
and the spin chirality, Therefore, the effective field becomes large when a light is emitted 
along z direction when a vortex structure is formed in the x-y plane. 

Next we study the magnetic vortex core switching by numerically solving the Landau- 
Lifshitz-Gilbert (LLG) equation including i/xiFE by using the forth-order Runge-Kutta 
method. As seen from Eq. ([5]), the effective field can not induce magnetization dynam- 
ics in zero temperature because the field direction is perfectly parallel (or anti-parallel) to 
the magnetization. The dynamics emerges when we include the thermal fluctuation of the 
magnetization, similarly to the case of the conventional IFE on a uniform ferromagnet^>^. 
The thermal fluctuation is included as a random field, C, whose ensemble average is defined 
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FIG. 2: (color online). Phase diagram of vortex core reversal for a disk with radius of 1 /xm and 
thickness of 20 nm at T = 100 K. L is proportional to the laser intensity. Region colored by blue 
represents negative value of averaged magnetization, (M^), indicating switching of a core, while 
red region represents no switching. 



where 7 is the gyromagnetic ratio and a is the Gilbert damping constant. We also include 
an internal field, needed to describe a vortex structure, representing the exchange in- 
teraction among the local spins and magnetic anisotropy. The LLG equation is therefore 
written as 



We consider a permalloy disk with a radius of i? = 1 ^m. We change the thickness of the 
disk from d = b nm to 40 nm. A vortex structure in the absence of light is calculated 
by minimizing the exchange energy and demagnetization energy. The radius of the vortex 
core is 25 nm. A circularly polarized light is irradiated perpendicular to the disk. The 
light-induced magnetic field, i^xiFE, is plotted in Fig. 1 (c). Even near a threshold value of 
L = 0.3 X 10^^^ Tm^, the maximum filed acting at the center of core reaches 0.3 T, but the 
filed decays rapidly as distance from the center increases. This field at the center is strong 
enough to induce picosecond core reversal. (In Ref. 2, uniform magnetic field of 0.03 T was 
used for a flip within 140 ps.) 



by a fluctuation-dissipation theorem as {Q)^ = and (C«(^)Ci(^'))c 




n =7 [n X {Htife + Hi + C)] 

- 07 [n X (n X (ifxiFE + Hi))] . 



(8) 
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FIG. 3: (color online). Thickness dependence of the switching time and the threshold intensity 
(L) for core switching. 

Figure 2 shows our result of averaged magnetization perpendicular to the vortex plane 
{i.e. the polarization of the vortex core) as a function of time and light intensity L = \L\. A 
circularly polarized light is irradiated at t = 0. We see immediately that there is a threshold 
value for the core flip at L ~ 0.3 x 10^^^ Tm^. (We defined the switching time time 
when the magnetization (M^) becomes —0.1 in Fig. 2 (see also Fig. 4 ).) The threshold 
value of L = 0.3 x 10^^^ Tm^ in Fig. 2 corresponds to the electric field of i? ^ 6 x 10^ V/m, 
and to the laser intensity of 10~^ J/m^ if a pulse duration is 0.25 ns. It is seen from Fig 2 
that a quite fast switching time (ts^ = 150 ps at L = 0.4 x 10"^^ Tm^) is realized. Compared 
with the current-induced vortex core switching, where rs„ = 20 ns the present mechanism 
is 100 times faster. Difference of the switching time between current-induced case and TIFE 
is understood by where there is a slow precessional motion. For the current-induce core 
switching process, a slow precession of a core is induced and a switching occurs when the 
core velocity becomes sufficiently large. On the other hand, TIFE does not induce such slow 
precessional motion. 

We found that the switching time t^^ depends much on the thickness d of the disk. As 
shown in Fig. 3, switching time becomes shorter when the disk becomes thinner. This is 
because the net magnetization to be reversed is proportional to the thickness. In contrast, 
the switching time is insensitive to the radius of the disk and pulse duration. 

From Fig. 2, one also notices that the switching time does not depend much on the 
applied laser power as long as it is above the threshold. In fact, even for much high intensity 
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FIG. 4: (Color online). Time-dependence of the averaged magnetization during the core switching 
process for several temperatures T = 0.1 K, 1 K, 10 K, and 100 K. The intensity of the light is 
L = 0.6 X 10~^^ Tm^. Inset: Dependence of switching time to temperature. 

of L = 6 X 10^^^ Tm~^, the switching time remains almost the same (~ 100 ps). This fact 
indicates that thermal fluctuation is essential to initiate the switching process, as noticed 
in a case of uniform ferromagnetr^. In fact, from the inset of Fig. 4, we see that the 
switching time depends on the temperature logarithmically. In the present TIFE, therefore, 
the polarized light intensity merely modifies the energy difference between the up and down 
states of vortex core, but does not affect the switching speed. This is because the magnetic 
field induced by the incident light, Eq. (jl]), is proportional to a local spin at each position. 
The field does not thus produce a torque which initiates the switching, and the initial torque 
must be induced by thermal fluctuation. 

In the conventional IFE in uniform ferromagnets, the switching time is inversely propor- 
tional to the temperature^. This dependence was explained by deriving an effective LLG 
equation by taking averaging of thermal field, C in Eq. ([8]). The effective Landau-Lifshitz- 
Gilbert equation for the average magnetic moment, m = (n), reads 

rh =7 [m X (iJiPE + -f^i)] - Tf^^n 

-a7([nx (nx (/fiFE + J/i))])c, (9) 

where i/jFE is the field generated by the conventional IFE and tl = ^^^^^ is a longitudinal 
relaxation time^. Since this relaxation time determines the switching time, the switching 
time of a uniform magnetization is inversely proportional to the temperature. In contrast 
to the conventional IFE on a uniform ferromagnet, TIFE for a vortex leads to a weak 
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logarithmic dependence on the temperature. This dependence would be understood by 
noting that spins in a vortex are strongly correlated to each other and thus a collective 
switching process not expressed by a single relaxation time occurs. 

In conclusion, we have discovered an unconventional inverse Faraday effect which appears 
when a circularly polarized light is irradiated to ferromagnet with a topological spin structure 
such as vortices and skyrmions. The effective magnetic field is derived analytically by 
using the non-equilibrium Green's function. We have shown by numerical simulation that 
a ultrafast switching of the magnetic vortex core is possible by use of TIFE resulting in a 
switching time of 150 ps, 100 times faster than that induced by the electric current. Even 
faster switching is expected if heating effect is included. Proposed mechanism opens a new 
way of the read-write process in spintronics devices. 
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